AMORTISSEMENT LANDAU 

C. MOUHOT AND C. VILLANI 

RÉSUMÉ. Dans cette note nous présentons les principaux résultats du récent tra- 
vail [14] , où le phénomène d'amortissement Landau est pour la première fois établi 
dans un contexte non linéaire. 

Abstract. In this note we présent the main results from the récent work [14] . 
which for the first time establish Landau damping in a nonlincar context. 

Mots-clés, amortissement Landau ; physique des plasmas ; astrophysique ; équation 
de Vlasov-Poisson. 

I. Introduction 

Le « modèle standard » de la physique des plasmas classique est l'équation de 
Vlasov-Poisson-Landau [HJEEE], ici écrite avec des conditions aux limites périodiques 
et en unités adimensionnées : 

(1) % + v- V x f + F[f] ■ VJ = ^ Q L (f, /), 

où f — f(t,x,v) est la fonction de distribution des électrons (t > 0, v G M 3 , 
x e T 3 = M 3 /Z 3 ), 



(2) F[f](t, x) = -jj ^W(x - y) fit, y, w) dwdy 

est la force auto- induite, W (x) = l/\x\ est le potentiel d'interaction coulomb ien, et 
Ql est l'opérateur de collision de Landau, décrit par exemple dans [TU] ou [T7] . Le 
paramètre A est très grand, variant typiquement entre 10 2 et 10 30 . 

Sur de très grandes échelles de temps (disons 0(A/ log A)), les phénomènes dissipa- 
tifs jouent un rôle non négligeable, et l'augmentation de l'entropie est supposée forcer 
la convergence (lente) vers une maxwellienne. Grâce aux progrès récents sur l'hy- 
pocoercivité, ce mécanisme est maintenant assez bien compris mathématiquement 
parlant, dès que l'on dispose d'estimations de régularité globale (voir [Ï6] et les 
références incluses). 
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Dix ans après avoir mis au point ce scénario collisionnel, Landau [9] formulait une 
prédiction beaucoup plus subtile : la stabilité d'un équilibre homogène vérifiant cer- 
taines conditions — par exemple fonction radiale de la vitesse v, mais pas forcément 
gaussienne — sur des échelles de temps beaucoup plus courtes (disons 0(1)), par le 
jeu de mécanismes purement conservatifs. Ce phénomène, appelé amortissement 
Landau, est une propriété de l'équation (non collisionnelle !) de Vlasov, obtenue en 
posant A = oo dans ([!]). La découverte de Landau fut jugée "stupéfiante" par ses 
contemporains ; c'est maintenant l'une des bases théorique de la physique classique 
des plasmas (parmi un très grand nombre de références citons [H [TU])- Des effets 
d'amortissement similaires ont été prédits en astrophysique [TT], mais aussi dans 
d'autres domaines de la physique. 

L'amortissement Landau est compris depuis longtemps au niveau de l'équation 
linéarisée [21 SI El El CES CE]. Cependant, il y a déjà cinquante ans de cela, dans le 
premier volume du Journal of Mathematical Physics, Backus [2] notait que l'échelle 
de temps de la linéarisation est parfois bien plus courte que l'échelle de temps de 
l'amortissement Landau, et exprimait des doutes sur la pertinence de l'approche par 
linéarisation. Son objection n'a jamais été réfutée, car les difficultées conceptuelles 
et techniques associées à l'équation non linéaire ont sévèrement limité les résultats 
obtenus dans ce cadre : dans [31 [S] on prouve seulement l'existence de certaines 
solutions amorties. 

Nous comblons cette lacune dans un travail récent [H] , dont nous décrivons main- 
tenant le résultat principal. 

2. RÉSULTAT PRINCIPAL 
Si / est une fonction définie sur T d x ~R d , on note, pour k G Z d , t] G M d , 

f(k,v)= ! f(x,v)e- 2l7Tkx alx, f(k,r))=(( f(x, v) e" 2 ^ e~ 2l ^ v alv dx. 

JT d J JT d xR d 

On pose en outre 

(3) \\f\\x,^ = sup(|/(M)l e 2 ^ e 2 ^) + [[ \f(x, v)\ e 2 ^ dvdx. 

k,ri v ' J JT^xR 1 * 

Théorème 1 (amortissement Landau non linéaire pour interactions générales). 
Soient d > 1, et f° : R — > M + un profil de vitesses analytique. Soit W : T d — ► K un 
potentiel d'interaction. Pour k G Z d , £ G C, on pose 

C(k,Ç) = -Air 2 W(k) / e 2 ^ H \f{kt)\\k\ 2 tdt. 

Jo 
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On suppose qu'il existe A > tel que 

~ \ n 

(4) sup \f°( V )\e 2 ^ < C , V -rWVîfWiAw < O), 

(5) inf inf \C(k,£) - il > « > 0, 

(6) 3 7 > 1; VA; G Z d ; \W(k)\ < 



dès que < A' < A, < // < \x, (3 > 0, r G N, il existe e > et C > 0, 
dépendant seulement de d, 7, A, A', /i, //, Cq, k, Cw, P, r tels que, si fi > vérifie 

(7) *:=||/i-/°||^<c, 
alors l'unique solution de l'équation de Vlasov non linéaire 

(8) ^ + vV x f + F[f]-V v f = 0, F[f](t,x) = - J VW(x-y)f(t,y,w)dwdy, 

définie pour tous les temps et telle que f(0,-) = fi, vérifie 

(9) \\p(t,-)-Poo\\ Cr{Td) <C5e- 2 ^\ 

où p(t,x) = f f(t,x,v)dv, = ff fi(x,v) dv dx. En outre il existe des profils 
analytiques f +OQ {y), f^^iy) tels que 

f(t, ■) t ^ ±QC > f ±oo faiblement 
f(t, x, ■ ) dx t ~* ± °° ) f ±oo fortement (dans C r (M^) ), 



ces convergences étant également en 0(ôe 27rA 'l'l). 

Ce théorème, entièrement constructif, est presque optimal comme le montrent les 
commentaires qui suivent. 

Commentaires sur les hypothèses : Les conditions aux limites périodiques 
sont bien sûr contestables; cependant, au vu des contre-exemples de Glassey et 
Schaeffer [5], un mécanisme de confinement — ou tout au moins de limitation des 
longueurs d'onde spatiales - - est indispensable. La condition (j3J) exprime quan- 
titativement l'analyticité du profil f°, sans laquelle on ne peut espérer avoir de 
convergence exponentielle. L'inégalité (jSJ) est une condition de stabilité linéaire, es- 
sentiellement optimale, qui couvre les cas physiquement intéressants : l'interaction 
de Newton aux longueurs d'onde inférieures à la longueur d'instabilité de Jeans, et 



4 C. MOUHOT AND C. VILLANI 

l'interaction coulombienne pour toutes les longueurs d'ondes autour de profils f° 
radialement symétriques en dimension supérieure ou égale à 3. La condition en 
revanche n'intervient que dans la stabilité non linéaire, elle inclue les interactions de 
Coulomb ou Newton dans le cas limite 7 = 1. Quant à la condition (J7|), son caractère 
perturbatif est naturel au vu des spéculations théoriques et des études numériques 
sur le sujet. 

Commentaires sur les conclusions 

(1) La convergence en temps grand est basée sur un mécanisme réversible, pu- 
rement déterministe, sans fonctionnelle de Lyapunov ni interprétation varia- 
tionnelle. Les profils asymptotiques /±oo gardent d'ailleurs la mémoire de la 
donnée initiale et de l'interaction. Cette convergence « sans raison d'être » 
n'était pas vraiment attendue, puisque la théorie quasilinéaire de l'amortis- 
sement Landau [H Vol. II, Section 9.1.2] ne prédit la convergence qu'après 
une prise de moyenne sur des ensembles statistiques. 

(2) On peut interpréter ce résultat dans l'esprit du théorème KAM : pour l'équation 
de Vlasov linéaire, la convergence est forcée par une infinité de lois de conser- 
vation qui rendent le modèle « complètement intégrable » ; dès que l'on ajoute 
un couplage non linéaire, les lois de conservation disparaissent mais la conver- 
gence demeure. 

(3) Étant donné un équilibre stable f°, un voisinage entier — en topologie ana- 
lytique — de f° est rempli par des trajectoires homoclines ou (en général) 
hétéroclines. C'est la dimension infinie qui permet ce comportement remar- 
quable de l'équation de Vlasov non linéaire. 

(4) La convergence en temps grand de la fonction de distribution n'a lieu qu'au 
sens faible ; les normes des dérivées en vitesse croissent très vite en temps 
grand, ce qui traduit une filamentation dans l'espace des phases, et un trans- 
fert d'énergie (ou d'information) des basses vers les hautes fréquences (« tur- 
bulence faible »). 

(5) C'est ce mécanisme de transfert d'information aux petites échelles qui permet 
de réconcilier la réversibilité de l'équation de Vlasov-Poisson avec l'apparente 
irréversibilité de l'amortissement Landau. Notons que le mécanisme « dual » 
de transfert d'énergie vers les grandes échelles, aussi appelé radiation, a été 
largement étudié dans le cadre des systèmes hamiltoniens de dimension infi- 
nie. 

On trouvera davantage de commentaires, aussi bien mathématiques que physiques, 
dans [H]. 
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3. Stabilité linéaire 

La stabilité linéaire est la première étape de notre étude ; elle ne demande qu'un 
investissement technique assez réduit. 

Linéarisée autour d'un équilibre homogène f°(v), l'équation de Vlasov devient 

ai p 

(10) — + v ■ V x h - (W * p) ■ VJ° = 0, p= hdv. 

Il est bien connu que cette équation se découple en une infinité d'équations indépendantes 
régissant l'évolution des modes de p : pour tout k G Z d et tout t > 0, 

(11) p{t,k)- f K Q {t-T,k)p{T,k)dT = hi{k,kt), 

Jo 

où hi est la donnée initiale, et K° un noyau intégral qui dépend de f° : 

(12) K°(t,k) = -A7r 2 W(k)f (kt) \k\ 2 t. 

On déduit alors de résultats classiques sur les équations de Volterra que pour tout 
k y£ la décroissance de p(t, k) quand t — > +oo est essentiellement contrôlée par le 
pire de deux taux de convergence : 

• le taux de convergence du terme source au membre de droite de (ITT1) . qui ne 
dépend que de la régularité de la donnée initiale dans la variable de vitesse ; 

• e~ xt , où À est le plus grand réel positif tel que la transformée de Fourier- 
Laplace (dans la variable t) de K° ne s'approche pas de 1 dans la bande complexe 
{0 < < A}. 

Le problème consiste donc à trouver des conditions suffisantes sur f° pour garantir 
la stricte positivité de À. Depuis Landau, cette étude est traditionnellement réalisée 
la formule d'inversion de la transformée de Laplace; cependant, dans la 
perspective de l'étude non linéaire, nous lui préférons une approche plus élémentaire 
et constructive, basée sur la simple formule d'inversion de Fourier. 

On montre ainsi l'amortissement Landau linéaire, sous les conditions dSJ) et dlj), 
pour n'importe quelle interaction W telle que VW G L 1 (T d ), et pour toute condition 
initiale analytique (sans restriction de taille dans ce contexte linéaire). On retrouve 
comme cas particulier tous les résultats précédemment connus sur l'amortissement 
Landau linéaire [H [T3J [15] ; mais on couvre également l'interaction newtonienne. En 
effet, la condition (jSJ) est vérifiée dès que l'une ou l'autre des conditions suivantes 
est vraie : 
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(a) VA; G Z d , M z G R, W(k) > 0, z<p' k (z) < 0, où <p k est la « marginale » de /° 
selon la direction k, définie par 

<P k (z)= [ f°(w)dw, zeR; 
J w\+ k± 

(b) 4vr 2 (max \W(k)\) \ sup f \f°(ra)\rdr) < 1. 

\\v\=iJo ) 

Cette dernière condition couvre l'interaction de Newton en-deçà de la longueur de 
Jeans. On renvoie à [TU Section 3] pour plus de détails. 

4. Stabilité non linéaire 

Pour traiter la stabilité non linéaire, nous commençons par introduire des normes 
analytiques « hybrides » (basées sur la taille des dérivées successives dans la va- 
riable de vitesse, et sur la taille des coefficients de Fourier dans la variable de position) 
et « glissantes » (la norme utilisée changera au cours du temps pour tenir compte 
des transferts dans l'espace des phases). Cinq indices permettent d'obtenir toute la 
souplesse nécessaire : 

(13) 11/11^), = EE^'^+w) 7 ^ 

k£Z d n£N d 

(Par défaut 7 = 0.) Un théorème fastidieux d'injection « à la Sobolev » compare ces 
normes à d'autres plus traditionnelles, telles que les normes ||/||a >( u,/3 de (j3J). 

Les normes Z possèdent des propriétés remarquables vis-à-vis de la composition 
et du produit. Le paramètre r permet de compenser, dans une certaine mesure, la 
filamentation. Enfin le caractère hybride est bien adapté à la géométrie du problème. 
Si / ne dépend que de x, la norme (fï3"l) coïncide avec la norme d'algbere jF Ar+/i ' 7 
définie par 

(14) = E e 2 ^ W (1 + \k\y. 

k&Z d 

(Nous utilisons également la version « homogène » jF Ar+/i ' 7 lorsque l'on n'inclut le 
mode k = dans la somme.) 

L'équation de Vlasov est résolue par un schéma de Newton dont la première 
étape est la solution du linéarisé autour de f° : 

f n = f + h 1 + . . . + h n , 



(V v + 2 ï 7TTk) n f(k, 



LP(dv) 
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n > 1, 



d t h n+1 
h n+1 (0, 



'dth 1 + v ■ V x h x + Flh 1 } ■ V v f° = 

fc 1 (o,-) = /<-/° 

v • V x ./i n+1 + F[/ n ] • V^ n+1 + F[/i n+1 ] ■ 
) = 0. 



V v f n = -F[h n ] ■ V v h n 



Dans un premier temps, on établit la régularité analytique en temps petit des 



h n (r, ■) en norme Z- 



A,(/*,7);i . 



; cette étape, dans l'esprit d'un théorème de Cauchy- 



Kowalevskaya, est réalisée grâce à l'identité 



(15) 



dt 



K 



1+T 



liv/|| 



■^(t).M(t)'iP) 



où X(t) =\-Kt, fj,(t) =n-Kt. 

Dans un deuxième temps, on établit des estimations uniformes en temps sur 
chaque h n , cette fois par une méthode partiellement eulérienne et partiellement la- 
grangienne, en intégrant l'équation le long des caractéristiques (X™ t , V™ t ) crées par la 
force F[f n ]. (Ici r est le temps initial, t le temps courant, (x, v) les conditions initiales, 
(X n ,V n ) l'état courant.) La régularité de ces caractéristiques est exprimée par des 
contrôles en norme hybride sur les opérateurs Çlf T (x, v) = (X™ T , Vt l T )(x + v(t — r),v), 
qui comparent la dynamique perturbée à la dynamique libre, informellement appelés 
opérateurs de scattering (en temps fini). 

On propage alors un certain nombre d'estimations à travers le schéma, les plus 
importantes étant (en simplifiant légèrement) 



(16) 
(17) 



sup 

T>0 



h n (r, ■ , v) dv 



sup 

t>r>0 



_A„(l+6),M»il — 

T 

1+6 



b(t) 



B 



? ^n(l + 6),(/Jn,7);°° 
1+6 



< c 



II 



6k e -2^(A fe -A„ +1 )T 



^ 2vr(A fc - A n+1 )^ 



min{t — r ; 1}. 



On note, dans ([16]) . l'augmentation linéaire de la régularité de la densité spatiale, 
qui fait contrepoids à la détérioration de régularité dans la variable de vitesse. Dans 
fTTTj) . le léger décalage des indices par la fonction b(t) sera crucial pour absorber les 
termes d'erreur provenant de la composition ; la constante B est elle-même choisie 
en fonction des estimations en temps petit réalisées précédemment. Enfin, dans (ÎT8l) . 
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remarquons le contrôle uniforme en t, et l'amélioration des estimées dans les deux 
régimes asymptotiques t — > r et r — > oo ; ceci également est important pour la 
gestion des termes d'erreur. Les constantes \ n et /i n décroissent à chaque étape du 
schéma, convergeant — pas trop rapidement — vers des limites Aoo, /i^ positives; 
dans le même temps, les constantes S n convergent extrêmement vite vers 0, ce qui 
garantit « par rétroaction » l'uniformité des constantes du membre de droite de (ÎT8l) . 

Les estimations f|T8|) sont obtenues par des applications répétées de théorèmes 
de point fixe en normes analytiques. Un autre ingrédient essentiel pour passer de 
l'étape n à l'étape n+1 est le mécanisme d'extorsion de régularité, que nous allons 
décrire dans une version simplifiée. Étant données deux fonctions de distribution / 
et /, dépendant de t, x, v, définissons 



a(t, x) 



f [ (F[f]-V v f)(r,x-v(t-r) } v)dvdr. 



L'interprétation de a est comme suit : si les particules distribuées selon / exercent 
une force sur des particules distribuées selon /, alors a est la variation de densité 
J f dv causée par la réaction de / sur /. Nous montrons que si / a une régularité 
glissante élevée, alors la régularité de a en temps grand est meilleure que ce que l'on 
attendrait : 

-t 



(19) 



ou 



K(t,r) 



\a(t, 



< 



[ K(t,r)\\F[f(r,-)] 
J o 



dr. 



Wvf(s, 



sup 

0<s<t 



1 + S 



1+T) SUp 



e -2n(\-\)\k(t-T)+£T\ e -2ir(p-n)\l\ 

l + \k-l\y 



Le noyau K(t,r) est d'intégrale 0{t) quand t — > oo, ce qui laisse craindre une in- 
stabilité violente ; mais au fur et à mesure que t augmente, il est également de plus 
en plus concentré sur des temps discrets r = kt/(k — t) ; c'est la manifestation des 
échos plasmas, découverts et observés expérimentalement dans les années soixante 
[12]. Le rôle stabilisant du phénomène d'écho, en relation avec l'amortissement Lan- 
dau, est mis à jour dans notre étude. Le mécanisme d'extorsion de régularité éclaire 
d'un jour nouveau l'interprétation qualitative populaire de l'amortissement Landau 
comme un transfert d'énergie d'onde vers particule. 

On analyse alors la réponse non linéaire due aux échos. Lorsque 7 > 11 on peut 
montrer que cette réponse est sous-exponentielle, de sorte que l'on peut la contrôler 
par une très légère perte de régularité glissante, au prix d'une constante gigantesque, 
qui sera plus tard absorbée par la convergence ultrarapide du schéma de Newton. 
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Finalement, une partie de la régularité de / aura été convertie en décroissance en 
temps long. 

Lorsque 7 = 1, des estimations plus fines sont nécessaires. Pour pouvoir traiter 
ce cas, nous étudions la réponse non linéaire mode par mode, c'est-à-dire en esti- 
mant f)(t, k) pour tous les k, via un système infini d'inégalités. Cela permet de tirer 
avantage du fait que les échoes qui se produisent à des fréquences différentes sont 
asymptotiquement plutôt bien séparés. Par exemple en dimension 1, l'écho dominant 
au temps t et à la fréquence k correspond à r = kt/(k + 1). 

En pratique, les trajectoires rectilignes dans (1191) doivent être remplacées par les 
caractéristiques (ceci traduit le fait que / exerce aussi une force sur /), ce qui est 
source de difficultés techniques considérables. Parmi les moyens mis en oeuvre pour 
les surmonter, mentionnons un deuxième mécanisme d'extorsion de régularité, en 
temps court et proche dans l'esprit des lemmes de moyenne ; en voici une version 
simplifiée : 

(20) \\a(t, •)\\j:\t+ lt < / \\F[f(r, ■ )] Il FA[T _ b(t _ r)]+M , 7 II V/(r, • ) Il a ( i +6 ), (m ,o);i dr. 

Jq ^ "^T-bt/il+b) 

On voit dans (J2Ô1 que la régularité de a est meilleure que celle de F[f], avec un gain 
qui dégénère en temps grand et lorsque r — > t. □ 
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